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Abstract 

The instability of protons is a crucial prediction of supersymmetric GUTs. We review 
the decay in minimal supersymmetric SU(5), which is dominated by dimension-five 
operators, and discuss the implications of the failure of Yukawa unification for the 



o 

. decay rate. In a consistent SU(5) model, where SU(5) relations among Yukawa cou- 

plings hold, the proton decay rate can be several orders of magnitude smaller than the 
present experimental bound. Finally, we discuss orbifold GUTs, where proton decay 
via dimension- five operators is absent. The branching ratios of dimension-six decay 
can significantly differ from those in four dimensions. 
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The stability of protons has been a research topic for a long time PP, and exactly 50 years 
ago, the first experiments started to measure its lifetime [21 — although no mechanism 
was known that leads to proton decay. The situation changed when both baryon number 
violating processes were found in the Standard Model (SM) jSj and the idea of Grand 
Unification came up, where the SM is embedded in a simple gauge group jl]. Today, it is 
widely believed that protons will decay, even if they do so after a tremendously long time. 

Grand Unified Theories predict proton decay and the unsuccessful search for decaying 
protons gives a strong constraint on GUT models. The first GUT model, minimal SU(5) 
could already be excluded in the early 1980s. Recently, it was claimed that even its 
supersymmetric version j3] is excluded by the experimental limit on proton decay [HI H] . In 
this review, we discuss these analyses and point out the crucial assumptions which have been 
made to analyze the dominant dimension-five operators. We will see that the problem of 
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minimal supersymmetric SU(5) is less proton decay than the failure of Yukawa unification. 
Hence, we need a consistent model which explains the Yukawa couplings and we will see 
that a supersymmetric SU(5) model with minimal particle content cannot be excluded by 
proton decay jH]. 

Although such models are in agreement with the experimental limits on proton decay, 
the dimension-five operators are troublesome, as the Yukawa couplings must obey special 
relations among themselves. One would therefore like to avoid these operators. Orbifold 
GUTs offer an attractive solution, since they forbid dimension-five operators and, moreover, 
enable us to deal with the doublet-triplet splitting problem P UHl HH 112] • Proton decay 
appears via dimension-six operators, which can lead to unusual final states depending on the 
localization of matter fields [THl ll4j. We will illustrate this on the basis of a six-dimensional 
SO(IO) model ^Hl- Hence, if proton decay is observed in the future, the measured branching 
ratios can make it possible to distinguish orbifold and four-dimensional GUTs. 



1 Minimal supersymmetric SU(5) 

We start this section by briefly describing the minimal supersymmetric SU(5) GUT model 
[S]. It contains three generations of chiral matter multiplets, 10 = (Q, m*^, 6°), 5* = {dP,L), 
and a vector multiplet A (24) which includes the twelve gauge bosons of the SM and twelve 
additional ones, the X and Y bosons. Because of their electric and color charges, the latter 
mediate proton decay via dimension-six operators. At the GUT scale, SU(5) is broken to 
GsM = SU(3) X SU(2) X U(1)y by an adjoint Higgs multiplet S(24). A pair of quintets, -f^(5) 
and if (5*), then breaks G^m to SU(3) x U(l)em at the electroweak scale. The superpotential 
is given by 

W = \miiT.^ + \aii + (S + 3ct) H + iF/^lOi ^^jH + V2Y^^\Qi 5*H (1) 
with m = O (Mgut)- The adjoint Higgs multiplet, 

acquires the vacuum expectation value (VEV) (S) = a diag (2, 2, 2, — 3, — 3) so that the X 
and Y bosons become massive, whereas the SM particles remain massless. The components 
Eg and S3 of S(24) both acquire the mass 

Ms = Mg = M3 = |m . (3) 

The pair of quintets, H{5) and H{5*), contains the SM Higgs doublets, Hf and Hj, which 
break GgM, and color triplets, and He, respectively. To have massless Higgs doublets 
Hf and Hf, while their color-triplet partners (leptoquarks) are kept super- heavy, the mass 
parameters of H{5) and H{5*) have to be fine-tuned O (^) ~ 10~^^. This is called the 
doublet-triplet-splitting problem. The RGE analysis gives constraints on the masses of the 
new particles [7]. 
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Figure 1: Proton Decay via dimension-five operators: They result from excliange of the leptoquarks 
followed by gaugino or higgsino dressing. 



Expressed in terms of SM superfields, the Yukawa interactions are 

Wy = Y^^ Hj + d!^ Hj + Y^^ L, Hf 

+ ly;^- Q, Ha + Yil L, + FJ^' «f + YZ H, 



(4) 



where 



— — Yue — ^1 , Y^ — Y^ — Y^i — Y^d — ^2 • (5) 



Apart from the SM couphngs, there are four additional ones due to the colored leptoquarks. 
Integrating out those leptoquarks, two dimension five operators remain which lead to proton 
decay (Fig. Hj) [E], 



M, 



He 



\ y;^ y'r (q.. q,) (Qk + («f e-) {ui 



(6) 



called the LLLL and RRRR operator. The scalars are transformed to their fermionic 
partners by exchange of a gauge or Higgs fermion. Neglecting external momenta, the triangle 
diagram factor reads, up to a coefficient k depending on the exchange particle, 

d^k I 1 1 1 

-/(M;mi,m2), (7) 



i(27r)4 ml-k'^ ml-k"^ M -li (47r) 



^/,^ X M / m? , m? mo , mi, i 



with 



ml - ml \m{ - ml - 

where M and denote the gaugino and sfermion masses, respectively. 

As a result of Bose statistics for superfields, the total anti-symmetry in the colour index 
requires that these operators are fiavor non-diagonal J7j. The dominant decay mode is 
therefore p — > Ki). Since the dressing with gluinos and neutralinos is fiavor diagonal, the 
chargino exchange diagrams are dominant fHl The wino exchange is related to the 
LLLL operator and the charged higgsino exchange to the RRRR operator, so that the 
coefficients of the triangle diagram factor are 

f^L = '2g'^ , Kj^ = yy' . (9) 
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Here y and y' denote the corresponding Yukawa couplings (cf. Fig. |l(b)| ) and g is the gauge 
couphng. 

Let us sketch how to calculate the decay rate via dimension-five operators; for details see 
Ref. jH]. The Wilson coefficients Csl = YqqYqi and C^^ = YueYud are evaluated at the GUT 
scale. Then they have to be evolved down to the SUSY breaking scale, leading to a short- 
distance renormalization factor Ag. Now the sparticles are integrated out, and the operators 
give rise to the effective four-fermion operators of dimension 6. The renormalization group 
procedure goes on to the scale of the proton mass, ~ 1 GeV, leading to a second, long- 
distance renormalization factor Ai [201 • At 1 GeV, the link to the hadronic level is made 
using the chiral Lagrangean method Thus the decay width can be written as 



/€ 1 

^had Ai -——f{M] mi, 1712) As — — Cs 
{Any Mh, 



(10) 



where the sum includes all possible diagrams. 

The decay width for the dominant channel p K^v reads 



2 2 \2 

mt — mi ^ 



, [ I'm — ^i^r I ^ — ^ 



3m R \ 3m fj 



+ Cf"" 1 - -^(3F-D' 
* 3mR 



(11) 



Here, mp and denote the masses of the proton and kaon, respectively, and /tt is the pion 
decay constant, ms is an average baryon mass according to contributions from diagrams 
with virtual S and A [21]. D and F are the symmetric and antisymmetric SU(3) reduced 
matrix elements for the axial- vector current. 

According to the two Wilson coefficients, the coefficients C5 split into two parts. 



C, = pC,^ + aC^^, (12) 



with 



Cll/rl = j^C,,/,,A,A^f{M;mi,m2) (13) 



and the hadron matrix elements a and /? 

oPlMp = ea/37 (0|(rfa^9 , pPj^Up = eai3y{0\{d^u'^)ul\p) . (14) 

While the hadronic parameters are fairly known, the masses and mixings of the SUSY- 
particles are unknown. We know through their absence that they have to be heavier than 
O (100 GeV); on the other hand, they are expected not to be much heavier than O (1 TeV). 
Looking at the dressing diagram we notice that when taking the sfermions to be degenerate 
at a TeV, the triangle diagram factor (jHJ is given by 

x/.^ ^ ^ f 2 1,^2 ,^2i M<t.m M 
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Figure 2: Decay rate T(p — > K~^i') as function of tan /? in minimal SU(5) with Ygi = Y^d = Y^- 




The simplest case is to assume that the sfermions have masses of 1 TeV. An exception is often 
made for top squarks. Since the off-diagonal entries of the mass matrix are proportional to 
m^, the mixing in the stop sector is expected to be large, with at least one eigenvalue much 
below 1 TeV. In analyses, one typically uses 400 GeV, 800 GeV, or 1 TeV for . For the other 
sfermions, the mixings are neglected. The proton decay rate is further suppressed by light 
gauginos and higgsinos. Note that the experimental limit for charginos is m^± > 67.7 GeV 

Since proton decay is dangerously large, the decoupling scenario ^24j has also been stud- 
ied, where the scalars of the first and second generation can be as heavy as 10 TeV [7] . Such 
an adjustment has been motivated by the supersymmetric flavor problem. In this scenario, 
proton decay via dimension-flve operators is clearly dominated by the third generation. 

Let us analyze the Wilson coefficients now. Minimal SU(5) predicts that the Yukawa 
couplings of down quarks and charged leptons are unified, as can be seen from Eqs. 
While mh = rrir can be fulfilled at the GUT scale, the equalities fail for the first and second 
generation. Nevertheless, proton decay via dimension-five operators has been analyzed 
assuming 

Yqq = Yue = Yu 5 Yql = Yud = Yd ■ (16) 

Then the decay rate can be calculated as described above. Note, however, that the choices 
Yqi = Yud = Yf, or Yqi = Yd, Y^d = Y^. would be equally justified. As we shall see, this 
ambiguity strongly affects the proton decay rate.^ 

Fig.|21shows the results of the following three cases: (i) all sfermions have masses of 1 TeV; 
(ii) mi is changed to 400 GeV; (iii) decoupling scenario, where the scalars of the first and 

^Another uncertainty concerns the sfermion mixings. Due to constraints by flavor changing neutral 
currents, they are assumed to coincide with the fermion mixings; see, however Refs. |25| . 
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Figure 3: Decay rate T{p K^v) as a function of tan P with Ygi = Y^d = Y^,. 



second generation have masses of 10 TeV. The dash-dotted hne represents the experimental 
hmit r = 6.7 X 10^^ years as given by the Super-Kamiokande experiment ISSIEE], the dotted 
hne is the newer hmit r = 1.9 x 10'^^ years [2Zj. 

The Wilson coefficients are proportional to tan [3 + ^^^^ and (tan (3 + ^j;^;^)^, where tan (3 
is the ratio of the vacuum expectation values of the Higgs doublets Hf and Hf. Clearly 
low values are preferred for obtaining a low decay rate; on the other hand, the top Yukawa 
coupling becomes non-perturbative for tan/? < 2.5. 

As already mentioned, it is possible to constrain the leptoquark mass using the RGEs. 
These constraints depend strongly on the Higgs representations, therefore we choose the 
most conservative value M^^ = Mqut = 2 x 10^^ GeV in order to study whether the 
experimental limit already rules out this model. 

The decay rate in Fig. |21is always above the experimental limit, which led to the claim 
that minimal SU(5) is excluded j^l E]. But as already discussed, there is no compehing 
reason for the assumption Ygi = Y^d = Y^l In order to illustrate the strong dependence of 
the decay rate on flavor mixing and therefore on Yukawa unification, let us study the case 

^gq ^ue ; ^ql ^ud ■ ("^''') 

Then the mixing matrix Ul Ud = Vckm, which appears in the Wilson coefficients, is replaced 
hj M. = UlUe- Note that the mixing matrix in Y^ or Yd is stih given by Vckm- Since Yd ^Yg, 
the masses and mixing of quarks and leptons are different and Ai is undetermined. 

We ffist ignore mixing, i.e. A4 = 1, and calculate the decay rate; the results are shown 
in Fig. 121 Without mixing, only scalars of the ffist and second generation take part, so the 
decay rate is significantly reduced in the decoupling scenario where the triangle diagram 
factor (jSI) changes by almost two orders of magnitude. 

If, moreover, we take Ai arbitrarily, it is possible to push the decay width below the 
experimental limit even for smaller sfermion masses. In the case = 400 GeV, this is only 
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possible for small values of tan/?. 

The fact that a sufficiently low decay rate can be found illustrates the dependence on 
flavor mixing and therefore the uncertainty due to the failure of Yukawa unification. Minimal 
supersymmtric SU(5) can only be excluded by the mismatch between the Yukawa couplings 
of down quarks and charged leptons as non-supersymmetric SU(5) is excluded by the failure 
of gauge unification. 



2 Consistent supersymmetric SU(5) 

Minimal supersymmetric SU(5) is inconsistent due to the failure of Yukawa unification, 
thus additional interactions are required which account for the difference of down quark 
and charged lepton masses. Such interactions are conveniently parameterized by higher 
dimensional operators, which are naturally expected as a result of interactions at a higher 
scale, where the GUT model is extended to a more fundamental theory. In particular, the 
GUT scale is only about two orders below the Planck scale. 

Including possible terms up to order 1/Mpi, the superpotential reads [2H1 

W = imtrS2 + iatrS=^ + 6^— ^ + c— — 
^ Mpi Mpi 

+ 7 eabcde f F/nof lOf i/^ + fl'lQflQf^Hf + /^nof lOf H''^ 
4 y ik/pi ■' ikZpi 

+ V2 (r^^HalOfS*, + K^a^lOf^b*, + h^H,10f^5*}j . (18) 

Now the masses of S3 and Sg are no longer identical, affecting the constraint on the lepto- 
quark mass PHI . 

The Yukawa couplings are then given by 

^ Mpi Mpi 

Here a/Mpi ~ O (10~^), and S and A denote the symmetric and antisymmetric parts of the 
matrices. Thus the three Yukawa matrices, which are related to masses and mixing angles 
at Mz by the RGEs, are determined by six matrices. From Eqn. (jl9p one reads off 

Yd-Ye = 5-^/i2 , (20) 

Mpi 

hence, the failure of Yukawa unification is naturally accounted for by the presence of h2. Note 
that no additional field is introduced to obtain this relation; it just arises from corrections 
0(a/Mpi). 
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Let us study the effect of those higher-dimensional operators on the Wilson coefficients. 
It is instructive to express the couplings in terms of the quark and charged lepton Yukawa 
couplings and the additional matrices / and /i, 

= y,', = y'e = - 5^ ( /f + ^/2") , = n + ^^h, , (21) 







u 




yA 


5 a 


u 


2 Mpi 



Note that Yqi — Yud = Yg — Yd, which means that Yqi and Yud cannot be zero at the same time. 
Nevertheless, proton decay via dimension-five operators can be avoided if both C^l = YqqYqi 
and C^R = YueYud vanish. For this purpose the couplings have to fulfill the relations 

/f + j/f = f^n'. /.'' = |^F/. (23) 

4 5 (7 5 cr 

This is only possible if we allow the (3,3)-component of /i and /2 to be O ('^) ^ 1. 
But even if we restrict ourselves to 'natural matrices', i.e. couplings up to we can 

considerably reduce the decay amplitudes. We will illustrate this with two simple examples 
where either the RRRR or the LLLL contribution vanishes at the GUT scale. 

We assume that Y^g, Yqi, Y^e and Yud are all diagonal by a suitable choice of matrices. 
The simplest form of Ygg and Y^e is then 

Ygg = Yue = dmgiO,0,yt) , (24) 

where yt is the top Yukawa coupling at Mqut- In the first model, we spread Y^ — Yd such 
that 

Yud = diag{0,ys-y^,,yb-yT) , Ygi = diag {ye - yd, 0,0) . (25) 

Clearly Cg-^™' = Y^i Y^^ is zero whenever a particle of the first generation takes part. Since 
at least one particle of the first generation is needed, the RRRR contribution vanishes com- 
pletely. After RGE evolution, this simple structure of Wilson coefficients changes slightly, 
but as shown in Fig. EJ the decay amplitude is always well below the experimental limit — 
in the case mi = 1 TeV even more than two orders of magnitude. Here, we restrict ourselves 
to the cases where the sfermion masses are not heavier than 1 TeV. 
If we choose the matrices Ygi and Y^d as 

Yud = diag{yd-ye,ys-y^^,yb) , Ygi = diag {0,0, yr) , (26) 

the LLLL contribution vanishes at M^ut because now 

^ijkm ^ ykrn ^^^j different 
from zero for i = j = k = 'm = 3, but the decay has to be non-diagonal. Only the RRRR 
contribution with a low absolute value remains. After renormalization, this contribution 
is still dominated by third generation scalars, and the LLLL operator contributes only via 
p —>■ Kut-. In this case, the decay rate is even smaller (Fig. EJ. Furthermore, due to the 
smaller (3,3)-component of hi compared to the first model, it can easily be used for higher 
values of tan /3. 
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Figure 4: Decay rate T(p K'^u) as function of tan /3 in the consistent models. 



We have seen that higher-dimensional operators can reduce the proton decay rate by 
several orders of magnitude and make it consistent with the experimental upper bound. This 
impressing fact leads to the question, whether there is any mechanism which would naturally 
lead to the required relations among Yukawa couplings. We can think of two possibilities, 
to start with some ad-hoc textures as a result of an unknown additional symmetry [22] or to 
extend the analysis to another group, in order to obtain additional symmetry restrictions. 
SO(IO), however, does not restrict the contributions from the higher-dimensional operators 
PU] . Thus it is more promising to consider theories, where the dimension-five operators are 
generically absent. 



3 Orbifold GUTs 

In orbifold GUTs, the GUT gauge symmetry is realized in more than four space-time di- 
mensions and broken to the standard model by compactification on an orbifold, utilizing 
boundary conditions that violate the GUT symmetry HDl HH 1121 ■ a consequence, the 
GUT and electroweak scale are separated in an elegant way and the dimension-five opera- 
tors are absent [101 E]- Furthermore, in case of SO(IO), the breaking to Ggm via the Higgs 
mechanism requires large Higgs representations and the path is not unique. Here, orbifold 
GUTs provide an attractive solution as well, as orbifold symmetry breaking can simplify 
the breaking pattern. 

We consider a 6D SUSY SO(IO) model compactified on a torus with three parities jHI]. 
The first Z2 breaks only = 1 6D SUSY to = 1 4D SUSY, the remaining ones SO(IO) 
as well, namely to Ggg = SU(5) x U(l) and Gps = SU(4) x SU(2) x SU(2) |S2], respectively. 
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Figure 5: The three SO(IO) subgroups at the corresponding fixed points (branes) of the orbifold 

Then the zero modes belong to the intersection of the two symmetric subgroups, GgM x U(l)'. 
The physical region is a 'pillow', corresponding to the two compact dimensions, with four 
fixed points as corners, where the unbroken subgroups are SO(IO), Gps, Gqg, and flipped 
SU(5) (Fig. ED- 

The field content is as follows [3^: The matter fields are located on the branes, whereas 
the gauge fields, six 10-dimensional Higgs representations and two combinations 16 © 16, are 
bulk fields. The matter fields include the 10 and 5* of SU(5) plus a singlet, the right-handed 
neutrino. With this setup, both the irreducible and reducible 6D gauge anomalies vanish. 

The main idea to generate fermion mass matrices is to locate the three sequential 16-plets 
on the three branes where SO(IO) is broken to its three GUT subgroups, namely ipi at Oqg, 
ip2 at Ofl and ip^ at Opg. The three 'families' are then separated by distances large compared 
to the cutoff scale = (9 (10^'' GeV), where the theory becomes strongly coupled. Thus 
they can only have diagonal Yukawa couplings with the bulk Higgs fields, direct mixings 
are exponentially suppressed. The brane fields mix with the bulk field zero modes. These 
mixings take place only among left-handed leptons and right-handed down-quarks, which 
leads to a characteristic pattern of mass matrices. 

The mass terms read 



W 



l„c_Ar c 
''3 



(27) 



Here and are diagonal 3x3 matrices but due to the mixing with the bulk field zero 
modes, m'', nf and are 4x4 matrices with a lopsided structure. 



tan (3 
1 

tan [3 



m ~ 



VuM^ 



diag (/ii,/i2,yU3) 
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III 










\ 
























/is 






V 


Ml 


M2 


M3 


M4 


/ 



(28) 



(29) 



fii and fij are O (v^w), whereas Ma are of the order of the compactification scale. 

The fermion masses and mixings agree well with the data within coefficients O (1) in the 
case /ii,2 — ^ [SI]. Note that while ~ m;,, the muon and electron masses can be easily 
accommodated because /ig and are not related by a fiipped SU(5) mass relation. 
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While the mixings for the left-handed down quarks and right-handed charged leptons 
are small so that their diagonalization matrices read 



(30) 



those for the right-handed down quarks and left-handed charged leptons are large. For 
those, the diagonalization matrices are given by 



( -Mx 

Ml2 

Ml 

Ml2 








p,M Mi2 

M2(/i3M3-/i3M4) 

fl M Mi2 



\ 



A3 Ml2 

M M 
A M 



Ml(A3M4+/i3M3) 
_ Af2(A3M4+M3M3) 



iM2 



Ml \ 

M 

M 
M 



A3 (m^+MI+MI ) -^^sMsMj 
uM2 



M4 
M 



(31) 



up to a two-dimensional mixing matrix for the second and third generation, which can be 
neglected for the purpose of proton decay. In Eqn. (p?T| . ft^ is (9 (/is, /is), furthermore 

Mf + M| and ~ 



Ml 
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The up-quarks are confined to one fixed point each, in particular the up quark is located 
on the Georgi-Glashow one. Therefore dimension-six proton decay can arise via the exchange 
of the SU(5) X and Y bosons as in the four- dimensional case. 



eflf 



9l 



^afi-y m Q fSi Q-yJ ^ ak'^ (3i Qyi -^fc] ~l~ h.C. . 



(32) 



The analysis of these operators is analogous to the dimension-five case; for details see 
Ref . [inj . Contrary to the 4D case, we have to take into account the presence of the Kaluza- 
Klein tower of the vector bosons V = {X, Y) with masses 



My{n,m) 



(2n + 1)2 ^ (2m)2 



Rl 



(33) 



Furthermore, the KK modes interact more strongly by a factor of -\/2 due to their bulk 
normalization [TTj . 

The sum over the KK modes is logarithmically divergent. Since the theory is only 
valid below the scale M^,, we restrict the sum to masses My(n,m) < M*. In the case 
i?5 = -Rg = 1/Mc, we obtain 



1 



TT 



^M2(ri,m) 



4M2 



+ 2.3 



(34) 



As we will see below, is constrained by the experimental limit of the dominant decay 
mode p — > e+7r° yielding M^/Mc — 12. 

Due to the parities and the SO(IO) breaking pattern, the coupling of the gauge bosons 
is not universal any more, in contrast to the 4D case. Their wavefunctions vanish on the 
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decay channel 


Branching Ratios [%] 




6D SO(IO) 


4D SU(5) 




case I 


case II 






75 


71 


54 




4 


5 


<1 


V 7r+ 


19 


23 


27 




1 


1 


<1 




<1 


<1 


18 




<1 


<1 


<1 



Table 1: Resulting branching ratios and comparision with SU(5) x \J{l)p. 



fixed points Opg and Ofj so that no coupling arises via the kinetic term with the second and 
third generation. The couplings to the bulk states are irrelevant since these are embedded 
in full SU(5) multiplets together with massive KK modes, thus the interaction always mixes 
the light states with the heavy ones jlSl HI] • 

Altogether, the operators for the decay via X and Y bosons read 



eS 



9l 



+ 2 e\ {pi^) vPa_\ dp„, [Ut] M71 e„/37 + h.c. (35) 



with the fermions in their mass eigenstates. 

We start with the simplest case of = and degenerate masses in the limit 
/^s = /^3) which we denote case I. In this case, the state (Iri has no strange-component and 
we obtain 



efl 



9l 



2 Vud "a M7 + J M7 e + 2 VudVus IJ''' Ua df3 



+ 2 VudVus e° ul + 2 V^^ /i^ < sp 

.7 = 1 



Vud d-y ~\~ Vus "^-y 



+ h.c. (36) 



The numerical results for the branching ratios are shown in Tabled 

Additional operators can arise at any brane from couplings which contain the derivative 
along the extra dimensions of the locally vanishing gauge bosons ^1]. Their contribution 
is suppressed with respect to the standard operators jJS]; the coefficients, however, are 
undetermined. For C5 = Cg = 1, the corrections for the listed branching ratios are less than 
3%. 

To make a comparison with ordinary 4D GUT models, we consider two SU(5) models 
described in Refs. jH3 tHH]; where we assume that some mechanism suppresses or avoids 
the proton decay arising from dimension-five operators.^ These models make use of the 

^For recent discussions of dimension-six operators in flipped SU(5) and SO(IO), see Refs. |37) . 
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Froggatt-Nielsen mechanism where a global U(l)^ symmetry is broken spontaneously 
by the vev of a gauge singlet field $ at a high scale. Then the Yukawa couplings arise from 
the non-renormalizable operators, 




Q 



Q. 



■ij ~ 9ij 



(37) 



Here, Qij are couplings O (1) and Qi are the charges of the various fermions. Particularly 
interesting is the case with a 'lopsided' family structure, where the chiral U(l)^ charges are 
different for the 5* and 10 of the same family. 

The difference between the 6D SO(IO) and 4D SU(5) models is most noticeable in the 
channel p ^'^K^. This is due to the absence of second and third generation weak eigen- 
states in the current-current operators and the vanishing (12)-component. Hence the decay 
p —>■ is doubly Cabibbo suppressed. This effect is a direct consequence of the local- 

ization of the 'first generation' to the Georgi-Glashow brane. 

Let us now consider the general case, where the M'^'^ are not degenerate and where /is 
and /ig'^ differ as well. From Eqn. (jHTj) we see that strange component in d!{ does not vanish 
anymore, but it is smaller than the bottom component. We have studied several cases 
whose results agree remarkably well. As an illustration, consider the case where /ig = 2/^3 
and /ig = 3/i3, with non-degenerate heavy masses Mf : M| : M| : Mf = \ '■ ^ 
Mf : M| : M| : Ml = | : ^ : 1 : | (case II). The branching ratios are listed in Tabled the 
differences between the two cases are indeed small. 

The most striking difference between the 4D and 6D case is the decay channel]? — >■ n~^K^, 
which is suppressed by about two orders of magnitude in the 6D model with respect to 4D 
models. It is therefore important to determine an upper limit for this channel in the 6D 
model. Varying the mass parameters in the range M^'^ / M = 0.1 — 1 and /ig''^/ /is = 0.1 — 10, 



one order of magnitude smaller than in the four-dimensional case. Note that in 5D orbifold 
GUTs, this channel is typically enhanced because the first generation is located in the bulk 



Finally, we can derive the lower limit on the compactification scale from the decay width 
of the dominant channel p e^n^ yielding Mc > 9 x 10^^ GeV, roughly half of the 4D GUT 
scale. 

4 Conclusion 

Supersymmetric GUTs provide a beautiful framework for theories beyond the Standard 
Model. The simplest model, minimal SU(5), is inconsistent due to the failure of down 
quark and charged lepton masses to unify. The decay amplitude therefore depends strongly 
on fiavor mixing. A consistent SU(5) model requires additional interactions, which are 



we find 



T{p e+vrO) 



<5%, 



(38) 



Ulllll. 



13 



conveniently parameterized by higher dimensional operators. These operators can reduce 
the proton decay rate by several orders of magnitude, hence proton decay does not rule out 
consistent supersymmetric SU(5) models. 

Nevertheless, it is interesting to consider models where dimension-five operators are 
absent; this is generically the case in orbifold GUTs. We discussed a 6D SO(IO) model, 
where the position of the matter fields on the different branes leads to a characteristic 
pattern of mass matrices. The branching ratios differ significantly from predictions of 4D 
GUTs, which can make it possible to distinguish orbifold and four dimensional GUTs. 
Furthermore, proton decay restricts the compactification scale Mc to be O (M^^^). 

If proton decay is observed in the future, there will be only a few events. The dominance 
of the channel p — > vK^ , however, would strongly point at dimension-five operators, whereas 
p e^TT^ refers to dimension-six decay. In the latter case, the absence of the process 
p fi'^K^ would disfavor the naive GUT model in four dimensions and test the idea that 
the different generations are spatially separated. 
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